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Abstract 



The so-called unified method expresses the solution of an initial-boundary value problem for 
an evolution PDE in the finite interval in terms of an integral in the complex Fourier (spectral) 
plane. Simple initial-boundary value problems, which will be referred to as problems of type I, 
can be solved via a classical transform pair. For example, the Dirichlet problem of the heat 
equation can be solved in terms of the transform pair associated with the Fourier sine series. 
Such transform pairs can be constructed via the spectral analysis of the associated spatial 
operator. For more complicated initial-boundary value problems, which will be referred to as 
problems of type II, there does not exist a classical transform pair and the solution cannot 
be expressed in terms of an infinite series. Here we pose and answer two related questions: 
first, does there exist a (non-classical) transform pair capable of solving a type II problem, 
and second, can this transform pair be constructed via spectral analysis? The answer to both 
of these questions is positive and this motivates the introduction of a novel class of spectral 
entities. Wc call these spectral entities augmented oigonfunctions, to distinguish them from 
the generalised eigenfunctions introduced in the sixties by Gel'fand and his co-authors. 



MSG: 35P10 (primary), 35C15, 35G16, 47A70 (secondary). 

1 Introduction 

Consider the following initial-boundary value problems for the linearized Korteweg-de Vries (LKdV) 

equation: 

Problem 1 



qt{x,t) + qxxx{x,t) 
q{x,0) 
q{0,t) = q{l,t) 

qxCi^t) 







9x(0,i)/2 







(x,f)G(o,i)x(o,r) 



t€[0,T], 
t€[0,T]. 



(1.1a) 
(1.1b) 
(1.1c) 
(l.ld) 



Problem 2 




{x,t) e {0,1) X {0,T), 
XG[0,1], 
t€[0,T]. 



(1.2a) 
(1.2b) 
(1.2c) 



q{x,0)^f{x) 
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Figure 1: Contours for the linearized KdV equation. 



It is shown in |3J [SJ [TOl [TT] that these problems are well-posed and that their solutions can be 
expressed in the form 

where, Fq is the circular contour of radius | centred at 0, F* are the boundaries of the domains 
{A e : Im(A^) > and |A| > 1} as shown on figure [l] a is the root of unity e^"^*/^, /(A) is the 
Fourier transform ^ 

/ e-'^''f{x)dx, AeC (1.4) 

and C^(^)j ^(-^) E^re defined as follows for all A e C: 
Problem 1 

A(A) = e'^ + ae'"^ + a^e'"'^ + 2(6"*^ + ae"*"^ + a^e"*"'^), (1.5a) 
(+(A) = /(A)(e'^ + 2ae-'"^ + 2a^e-"^^^) + f{a\){ae"^^ - 2ae-'^) 

+ /VA)(a'e*"'^ - 2a"e-'^), (1.5b) 
C-(A) = -/(A)(2 + a^e-"'^ + ae"^"'^) - a/(aA)(2 - e"'"'^) - a^f (a'X){2 - e"'"^). (1.5c) 
Problem 2 

A(A) = e-'^ + ae-'°-^ + o?er''^^^, (1.6a) 

C+(A) = /(A)(ae-*"^ + a^e"*"'^) - (a/(aA) + aVVA))e-'^, (1.6b) 

r(A) = -/(A) - a/(aA) - a2/(a'A). (1.6c) 

For evolution PDEs defined in the finite interval, x G [0, 1], one may expect that the solution 
can be expressed in terms of an infinite series. However, it is shown in [51 llOj that for generic 
boundary conditions this is impossible. The solution can be expressed in the form of an infinite 
series only for a particular class of boundary value problems; this class is characterised explicitly 
in |10j . In particular, problem 2 does not belong to this class, in contrast to problem 1 for which 
there exists the following alternative representation: 

E I '•^-^""Sf + ^ E / ."<.-««...r(Ai,,_ ,,,, 

A(a)=0 A(cr)=0 
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where To- is a circular contour centred at a with radius \\ the asymptotic formula for a is given 
in By using the residue theorem, it is possible to express the right hand side of equation (1.7) 
in terms of an infinite series over <j. 

We note that even for problems for which there does exist a series representation (like prob- 
lem 1), the integral representation (1.3) has certain advantages. In particular, it provides an 
efficient numerical evaluation of the solution [3] . 

Generic initial-boundary value problems for which there does not exist an infinite series repre- 
sentation will be referred to as problems of type II, in contrast to those problems whose solutions 
possess both an integral and a series representation, which will be referred to as problems of type I, 

existence of a series representation: type I 
existence of only an integral representation: type II. 



Transform pair 

Simple initial-boundary value problems for linear evolution PDEs can be solved via an appropriate 
transfrom pair. For example, the Dirichlet and Neumann problems of the heat equation on the 
finite interval can be solved with the transform pair associated with the Fourier-sine and the 
Fourier-cosine series, respectively. Similarly, the series that can be constructed using the residue 
calculations of the right hand side of equation (1.7) can be obtained directly via a classical transform 
pair, which in turn can be constructed via standard spectral analysis. 

It turns out that the unified method provides an algorithmic way for constructing a transform 
pair tailored for a given initial-boundary value problem. For example, the integral representa- 
tion (1.3) gives rise to the following transform pair tailored for solving problems 1 and 2: 



/(x) ^ F{\) : 
F(A) ^ fix) : 



'j^q^+{x,X)fix)dx ifAer+UFo, 

J^q)-ix,\)f{x)dx ifAer-, 



r+ 



r- 



(1.8a) 
(1.8b) 



where for problems 1 and 2 respectively, 4>^ are given by 
1 



(/.+ (a;. A) 



0'(a;, A) 



and 



0+(x,A) 
(l)^{x,X) 



27rA(A) 



27rA(A) 



1 

27rA(A) 
27rA(A) 



-iXx { AX 



2ae 



-iaX 



2a"e-*" ^) -f e-'"^^(ae 

I ^ — ia'^Xx/^2Aa'^X 



iaX 



-iXx 



(2 -t- a^e-'"-' -f 



ae 



-i\x 



{ae 



-iaX 



■ a e 



-ia-'x 



) - (ae-^^-^^^ + 



a e 



-ia Xx 



-iX 



(1.8c) 

(1.8d) 

(1.8e) 
(1.8f) 



The alternative representation (1.7) gives rise to the following alternative transform pair tailored 
for solving problem 1: 



F(A) ^ fix) 



(1.9) 



erg 
A(ct)=0 
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where -Fa(/) is defined by equations (1.8a), (1.8c) and (1.8d) and is defined below (1.7) 



The vaUdity of these transform pairs is estabhshed in section [2] The solution of problems 1 
and 2 is then given by 

qix,t)^U(^e^^''F^if)). (1.10) 

Spectral representation 

The basis for the classical transform pairs used to solve initial-boundary value problems for linear 
evolution PDEs is the expansion of the initial datum in terms of appropriate eigenfunctions of the 
spatial differential operator. The transform pair diagonalises the associated differential operator in 
the sense of the classical spectral theorem. The main goal of this paper is to show that the unified 



method yields an integral representation, like (1.3 1, which in turn gives rise to a transform pair 



like (1.8), and furthermore the elucidation of the spectral meaning of such new transform pairs 
leads to new results in spectral theory. 

In connection with this, we recall that Gel'fand and coauthors have introduced the concept of 
generalised eigenfunctions [B] and have used these eigenfunctions to construct the spectral repre- 
sentations of self-adjoint differential operators [51. This concept is inadequate for analysing the 
IBVPs studied here because our problems are in general non-self-adjoint. Although the given 
formal differential operator is self-adjoint, the boundary conditions are in general not self-adjoint. 

In what follows, we introduce the notion of augmented eigenfunctions. Actually, in order to 
analyse type I and type II IBVPs, we introduce two types of augmented eigenfunctions. Type I 
are a slight generalisation of the eigenfunctions introduced by Gel'fand and Vilenkin and are also 
related with the notion of pseudospectra [51. However, it appears that type II eigenfunctions 
comprise a new class of spectral functionals. 

Definition 1.1. Let C be a linear topological space with subspace $ and let L : ^ ^ C be a linear 
operator. Let 7 be an oriented contour in C and let E = {Ex : A G 7} be a family of functionals 
E\ G C". Suppose there exist corresponding remainder functionals Rx G and eigenvalues 
z : 7 — >■ C such that 

Ex{Lcp) = z{\)Ex{4>) + Rx{4>), V0G$, VAG7- (l-H) 



/"e'^^i?A(0)dA = O, V(/)G VxG [0,1], (1.12) 
J7 



then we say E is a family of type I augmented eigenfunctions of L up to integration along 7. 



i?^(0)dA = O, V G $, V a; G (0, 1), (1.13) 
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z{\y 

then we say E is a family of type II augmented eigenfunctions of L up to integration along 7. 



We note that the class of families of augmented eigenfunctions of a given operator is closed 
under union. 

Recall that in the theory of pseudospectra it is required that the norm of the functional Rx{<P) is 
finite, whereas in our definition it is required that the integral of eyi\){i\x)Rx{4') along the contour 
7 vanishes. Recall that the inverse transform of the relevant transform pair is defined in terms of 
a contour integral, thus the above definition is sufficient for our needs. 

It will be shown in Section [4] that {Fx : A G 3 cr G C : A((t) = 0} is a family of type I 
augmented eigenfunctions of the differential operator representing the spatial part of problem 1 
with eigenvalue A"^. Similarly {Fx : A G Fq} is a family of type I augmented eigenfunctions of the 



spatial operator in problem (1.2). However, {Fx : A G r+ U F } is a family of type II augmented 
eigenfunctions. 



4 



Diagonalisation of the operator 

Our definition of augmented eigenfunctions, in contrast to tlie generalized eigenf unctions of Gel'fand 
and Vilenkin Section 1.4.5], allows the occurence of remainder functionals. However, the con- 
tribution of these remainder functionals is eliminated by integrating over 7. Hence, integrating 
equation ( |1.11[ ) over 7 gives rise to a non-self- adjoint analogue of the spectral representation of an 
operator. 

Definition 1.2. We say that E = {Ex : A G 7} is a complete family of functionals Ex e C if 

(pe<^> and Ex(f> = Oy XE'^ =i> = 0. (1.14) 

We now define a spectral representation of the non-self- adjoint differential operators we study 
in this paper. 

Definition 1.3. Suppose that E — {Ex : A G 7} is o system of type I augmented eigenfunctions 
of L up to integration over 7, and that 

1 e'^'^ExLcjydX converges V V a; e (0, 1). (1.15) 

Furthermore, assume that E is a complete system. Then E provides a spectral representation of 
L in the sense that 

( e'^^'ExLcfxiX^ I e'^''z{\)Ex4>d\ V V a; G (0, 1). (1.16) 

Definition 1.4. Suppose thatE^^^ — {Ex '■ A G 7^^-'} is a system of type I augmented eigenfunctions 
of L up to integration over 7'^^ and that 

[ e'^^'ExLcjidX converges V € $, V a; G (0, 1). (1.17) 

Suppose also that ^ {Ex : X e 7(")} is a system of type II augmented eigenfunctions of L up 
to integration over j^^^"^ and that 

I e'^^Ex(l)dX converges V G $, V a; e (0, 1). (1.18) 

Furthermore, assume that E = u is a complete system. Then E provides a spectral 
representation of L in the sense that 

I e'^''ExL(j)dX= [ z{X)e'^''Ex(j)dX V € V a; € (0, 1), (1.19a) 

/ -^e'^''ExL(j)dX= [ e'^'=Ex(l)dX V G V a; e (0, 1). (1.19b) 

7^(11) z(Xj J^(ii) 

According to Definition 1.3 the operator L is diagonalised (in the traditional sense) by the 
complete transform pair 

Ex, f e'^'' ■ dx] . (1.20) 
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Hence, augmented eigenfunctions of type I provide a natural extension of the generalised eigen- 
functions of Gel'fand & Vilenkin. This form of spectral representation is sufficient to describe the 
transform pair associated with problem 1. However, the spectral interpretation of the transform 
pair used to solve problem 2 gives rise to augmented eigenfunctions of type II, which are clearly 
quite different from the generalised eigenfunctions of Gel'fand & Vilenkin. 

Definition 1.4 describes how an operator may be written as the sum of two parts, one of which 
is diagonalised in the traditional sense, whereas the other possesses a diagonalised inverse. 
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Theorem 1.5. The transform pairs {F\, fx) defined in (1.8a|-(1.8d) and (1.8a), (1.8b|, (1.8e) 



and (1.8f| provide spectral representations of the spatial differential operators associated with prob- 
lems 1 and 2 respectively in the sense of Definition\1.4\ 



Theorem 1.6. The transform pair {F\, f]^) defined in (1.8a), (1.8c), (1.8d) and (1.9) provides a 
spectral representation of the spatial differential operator associated with problem 1 in the sense of 
Definition | j.^t 

Remark 1. Both problems 1 and 2 involve homogeneous boundary conditions. It is straight- 
forward to extend the above analysis for problems with inhomogeneous boundary conditions, see 
Section [3. II 



2 Validity of transform pairs 



In section 
section 12.2 



2.1| we will derive the validity of the transform pairs defined by equations (1.8). In 
we derive an analogous transform pair for a general IBVP. 



2.1 Linearized KdV 



Proposition 2.1. Let F\{f) and fx{F) be given by equations (1.8a)-(1.8d). For all f <E C° 

such that /(O) /(I) ^ and /'(O) 2/'(l) and for all x G (0, 1), we have 

fAFxif)) = fix). 



[0,1] 



(2.1) 



Let Fxif) and /^(F) be given by equations ( [TSal , pTSbf , ( [L8c| and ( (LS^ . For all f e C°°[0,1] 
such that /(O) = /(I) = /'(I) = and for all x G (0, 1), 



fAFxif)) = fix). 



Proof. The definition of the transform pair ( 1.8a)-( 1.8d) implies 



fAFxif)) 



1 

2tt 



' AiX) 



1 

2tt 



i\{x-l) C (^) 

A(A) 



dA, 



(2.2) 



(2.3) 



where and A are given by equations ( 1.5 ) and the contours , T~ and Fq are shown in figure[T] 
The fastest-growing exponentials in the sectors exterior to are indicated on figure [2^. Each 
of these exponentials occurs in A and integration by parts shows that the fastest-growing-terms 
in are the exponentials shown on figure [2^ multiplied by X~^. Hence the ratio C^(A)/A(A) 
decays for large A within the sector tt/3 ^ argA ^ 27r/3 and the ratio C^(A)/A(A) decays for 
large A within the sectors —n ^ argA ^ — 27r/3, —tt/3 ^ argA ^ 0. The relevant integrands 
are meromorphic functions with poles only at the zeros of A. The distribution theory of zeros of 
exponential polynomials |7] implies that the only poles occur within the sets bounded by F^. 

The above observations and Jordan's lemma allow us to deform the relevant contours to the 
contour 7 shown on figure[2]D; the red arrows on figure[2^ indicate the deformation direction. Hence 
equation (2.3) simphfies to 



fAFxif)) 



1 

2^ 



A(A) 



(C+(A)-e-^\-(A)) dA 



Equations (1.5) imply, 



(C+(A)-e-\-(A))=/(A)A(A), 



(2.4) 



(2.5) 



where / is the Fourier transform of a pieccwisc smooth function supported on [0,1]. Hence the 
integrand on the right hand side of equation (2.4) is an entire function, so we can deform the 



contour onto the real axis. The usual Fourier inversion theorem completes the proof. 



The proof for the transform pair (1.8a), (1.8b), (1.8e) and (1.8f) is similar 



□ 
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Although 



Figure 2a Figure 2b 

Figure 2: Contour deformation for the hnearized KdV equation. 



(2.6) 



the values at the endpoints can be recovered by taking apporpriate limits in the interior of the 
interval. 



2.2 General 

Spatial diflferential operator 

Let C — C°°[0, 1] and : C — > C be the following linearly independent boundary forms 

n-l 

Bj<ty = J2 {bjkcl^^'Ho) + Pjk<l>^'Hl)) , J e {1,2, . . . 

A;=0 



(2.7) 



Let $ {0 e C : Bjcf) = V j G {1, 2, . . . , n}} and {B* : j e {1,2,..., n}} be a set of adjoint 
boundary forms with adjoint boundary coefficients 6*^,, Let 5 : $ — > C be the differential 

operator defined by 

d" 



(2., 



Then S is formally self-adjoint but, in general, does not admit a self-adjoint extension because, in 
general, Bj ^ B* . Indeed, adopting the notation 

n-l 

im^) = Hr ^(-l)^'(<^("-^--''Ha;)V;'-''\x)), (2.9) 
of [U Section 11.1] and using integration by parts, we find 

((-zd/dx)"0,^) = [cj^m) - + i-id/dxr^p), v 0, e c^io, i]. (2.10) 

If g $, then ^/j must satisfy the adjoint boundary conditions in order for [^V-'Kl) ~ [0'0](O) — 
to be valid. 
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Initial- boundary value problem 

Associated with S and constant a G C, we define tiie following homogeneous initial-boundary value 
problem: 

{dt + aS)q{x,t) =0 V (a;,0 G (0,1) X (0,r), (2.11a) 

g(x,0)=/(x) Va;e[0,l], (2.11b) 

9(-,i)e* Vt€[0,l], (2.11c) 

where / G $ is arbitrary. Only certain values of a are permissible. Clearly a = is nonsensical 
and a reparametrisation ensures there is no loss of generality in assuming \a\ — 1. The problem 
is guaranteed to be ill-posed (for the same reason as the reverse-time heat equation is ill-poscd) 
without the following further restrictions on a: if n is odd then a ~ ±i and if n is even then 
Re(a) ^ mHO]. 



A full characterisation of well-posedness for all problems (2.11) is given in |8l[T01[TT]; For even- 
order problems, well-posedess depends upon the boundary conditions only, but for odd-order it is 
often the case that a problem is well-posed for a = i and ill-posed for a = — i or vice versa. Both 



problems (1.1 1 and (1.2 1 are well-posed 



Definition 2.2. We classify the IBVP (2.11) into three classes using the definitions of 11 If : 



type I: if the problem for {S, a) is well-posed and the problem for (S, —a) is well-conditioned. 
type II: if the problem for {S, a) is well-posed but the problem for (5, —a) is ill-conditioned. 
ill-posed otherwise. 

We will refer to the operators S associated with cases I and II as operators of type I and type II 
respectively. 

The spectral theory of type I operators is well understood in terms of an infinite series repre- 
sentation. Here, we provide an alternative spectral representation of the type I operators and also 
provide a suitable spectral representation of the type II operators. 

Transform pair 

Let a = e^'^*/". We define the entries of the matrices Af^(A) entrywise by 



Mk,W - 5] (-*«'=-! A)'-&*,, (2.12a) 

r=0 
n-1 

Ml.^{\) = Y.{-ia^-^\rP],. (2.12b) 

r=0 

Then the matrix Af(A), defined by 

M,.,(A) = (A) + M-.(A)e— (2.13) 

is a realization of Birkhoff 's adjoint characteristic matrix. 

We define A(A) = detA/(A). From the theory of exponential polynomials [7], we know that 
the only zeros of A are of finite order and are isolated with positive infimal separation 5e, say. We 
define X^^ as the [n — 1) x (n — 1) submatrix of M with (1, 1) entry the (/ + 1, j + 1) entry of M. 
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The transform pair is given by 



F+if) ifAer+ur+ 



f{x)^F{X): Fx{f) = < i V (2-14a) 

[F^if) if Aer+ur^, 



^(A) ^ fix) : /.(F) 
where, for A e C such that A(A) 7^ 0, 



J e'^''F{X)dX, a;€[0,l], (2.14b) 



n n „ 1 

^^detX'^(A)Mf^.(A) / e-'°'"'^V(x)dx, (2.15b) 
;=i ,=1 -^0 



27rA(A) 

and the various contours are defined by 



r - To U r„ (2.16a) 

ro = r+ur~, (2.i6b) 

ro+= U ^(^'^)' (2.16c) 



A(a)=0 



^0= U ^(^'^)' (2.16d) 



A(cr)=0 



r„-r+ur-, (2.i6e) 

is the boundary of the domain 

{A e C± : Re(aA") > 0} \ |J D{<t, 2e). (2.16f) 



A(ct)=0 



Figure [3] shows the position of the contours for some hypothetical A with zeros at the black 
dots. The contour Tq is shown in blue and the contour Fq is shown in black. The contours Fj 
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and r~ are shown in red and green respectively. This case corresponds to a ~ —i. The figure 
indicates the possibility that there may be infinitely many zeros lying in the interior of the sectors 
bounded by F^. For such a zero, F^ has a circular component enclosing this zero with radius 2e. 
The validity of the transform pairs is expressed in the following proposition: 

Proposition 2.3. Let S be a type I or type II operator. Then for all f £ and for all x £ (0, 1), 



UFxif)) = 




e^^-F+{f)dX + 



Proof. A simple calculation yields 




e^^^F^{f)dX = f{x). 



(2.17) 



V/eC, V5, F+(/)-F-(/) 



2n 



/(A). 



(2.18) 



As shown in [TD], the well-posedness of the initial-boundary value problem implies F^{f) — 
0{X~^) as A — >■ oo within the sectors exterior to F^. The only singularities of F^{f) are isolated 
poles hence, by Jordan's Lemma, 




cret 
Iin(cr)>£, 
A(ct)=0 



iXx 77' + 




C(cr,£) JC(a,2e) 



'FUf)dX+ J2 

Im(cr)<£, 
A(a)-0 




Fxif)dX 



C{a,2€) 



{F+{f)~F^{f)) dX, (2.19) 



where 7 is a contour running along the real line in the increasing direction but perturbed along 
circular arcs in such a way that it is always at least s away from each pole of A. The series on 



the right hand side of equation (2.19) yield a zero contribution. As / € its Fourier transform / 



is an entire function hence, by statement (2.18), the integrand in the final term on the right hand 



side of equation (2.19) is an entire function and we may deform 7 onto the real line. The validity 
of the usual Fourier transform completes the proof. □ 



3 True integral transform method for IBVP 



In section 3.1 we will prove equation ( 1.10 ) for the transform pairs ( 1.8 ). In section 3.2 we establish 



equivalent results for general type I and type II initial-boundary value problems. 



3.1 Linearized KdV 



Proposition 3.1. The solution of problem 1 is given by equation (1.10), with F\{f) and fx{F) 
defined by equations (1.8a)-(1.8d). 



The solution of problem 2 is given by equation (1.10), with Fx{f ) and fx{F) defined by equa- 
tions (Oal, (Obi), pTSel) and (Ol. 



Proof. We present the proof for problem 2. The proof for problem 1 is very similar. 

Suppose q e C°°([0,1] X [0,T]) is a solution of the problem (1.2). Applying the forward 
transform to q yields 



Fx{q{;t)) 



_Ul(l)+{x,X)q{x,t)dx ifAeC+, 



{x^X)q{x^t) dx if A e 



(3.1) 



10 



The PDE and integration by parts imply the following: 



dt 



Jo 

-dlq{l, t)(t>^il, A) + dlq{0, i)(/)±(0, A) + a,(7(l, (1, A) - 9,^(0, t)5,(/.±(0, A) 

- q{l, t)d,,^^{l, A) + g(0, i)a,,</.±(0, A) + iX^'Fxiqi; t)). (3.2) 



Rearranging, multiplying by e '"^ * and integrating, we find 

2 

Fx{q{;t)) = e*^'*i^;,(/) + e*^'*X](-l)^- [^^^.^^(O, A)g,(0, A) - S^^<^±(1, A)0,(1, A)] , (3.3) 

3=0 



where 



ix,X)= / e-'^ 'diqix,s)ds. 



Evaluating 9^^='=(0, A) and 9^^^(1,A), we obtain 



Fx{q{.,t)) = e^''Fx{f) + 



21T 

+Qo(0,A)A2 
+(3o(l,A)A 



gi(l,A)iA(a-a2) 



A(A) 



A(A) 

2 (1 - a2)e'"^ + (1 - Q;)e*«'-^ 



A(A) 



+ (32(0,A) + Qi(0,A)iA 



for all A e C+ and 



(3.4) 



(3.5) 



Fx{q{.,t)) = e'^''Fx{f) + ^ 



iX't 



27r 



gi(l,A)iA 



+go(0,A)A' 



A(A) 



^ ^ go(i,A)A^-'^ + -''^^ + ^ 



A(A) 



A(A) 



+ g2(l,A)e- 



-iX 



(3.6) 



for all A e C". 
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Hence, the validity of the transform pair (Proposition 2.1 1 imphes 

iXx+rX'tp^^j^ dA 



qix,t) = <j / + / + 

Jr+ Jr 



?o(0,A)A- 



1 

+ 2^ 



Qi{l,X)iX{a-a^) 



A(A) 



1 

2t: 



A(A) 

iXx+iX^t 



?o(l,A)A 



r- 



2(l-a2)e'"-^ + (l-a)e*" ^ 

ia^X 



dA 



+Qo(0,A)A2 



A(A) 



gl(l,A)^A 
Qo(1,A)A2 



A(A) 



1 

2tt 



r+ 



A(A) 

^^-+^^'*[Q2(0,A) + Qi(0,A)iA] dA 



dA 



1 

2n 



by r±. Hence, 
equation (3.7) vanish. The boundary conditions imply 

Qo(0,A) = go(l,A) = Qi(l,A) -0, 



e'^("-i)+*^'*Q2(l,A)dA. (3.7) 

(3.8) 



Integration by parts yields 

as A — > oo within the region enclosed by r±. Hence, by Jordan's lemma, the final two lines of 

(3.9) 



so the second, third, fourth and fifth lines of equation (3.7) vanish. Hence 

f / + / |e*^"+'^'*FA(/)dA. 



q{x,t) 



(3.10) 

□ 



The above proof also demonstrates how the transform pair may be used to solve a problem 
with inhomogeneous boundary conditions: consider the problem, 



qt{x,t) + qxxxix,t) = 

q{x,0) = 0(a;) 
q{0,t) = hi{t) 
q{l,t)^h2{t) 
fc(l,i) = h^it) 



(x,t) e (0,1) X (0,T), 
a^e [0,1], 

[o,r], 
[o,r], 
[o,r]. 



(3.11a) 
(3.11b) 
(3.11c) 
(3.11d) 
(3. lie) 



for some given boundary data hj e C°°[0, 1]. Then Qo{Q, A), Qo{l, A) and Qi(l, A) are nonzero, but 
they are known quantities, namely i-transforms of the boundary data. Substituting these values 
into equation (3.7) yields an explicit expression for the solution. 

3.2 General 

Proposition 3.2. The solution of a type I or type II initial-boundary value problem is given by 



qix,t)^f,[e-^^"'F^if)). 



(3.12) 
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Lemma 3.3. Let / G $. Then there exist polynomials Pj of degree at most n — 1 such that 



F+(5/) = A"F+(/) + P+(A), 
F^iSf)^X"F^{f) + Py{X)e-^\ 



(3.13a) 
(3.13b) 



Proof. Let (0, V') be the usual inner product (j){x)^{x) dx. For any A € F, we can represent 
as the inner product F^{f) = (/, (/)J), for the function ^^(a;), smooth in x and meromorphic in A, 
defined by 



4>x {x) 



27rA(A) 

^ ' 1=1 .1=1 



^^detX'^"(A)Af+(A)e— 



27rA(A) 

^ ' 1=1 .1=1 



^^detX'^(A)Afr,(A)e-»'"'^^ 



As (j)^, Sf e C°°[0, 1] and a^'-^)" = 1, equation ( plo] ) yields 

^^±(5/) ^ A"F±(/) + [/0±](1) - [/0±](O). 
If B, B* : C°°[0, 1] C", are the real vector boundary forms 

B = {Bi,B2, . . . , Bn), B* = {B*, i?2, . ■ . , B*) 



(3.14a) 
(3.14b) 

(3.15) 
(3.16) 



then the boundary form formula [1, Theorem 11.2.1] guarantees the existance of complimentary 
vector boundary forms B^ B* such that 



[f^^m - [f^tm = Bf ■ B:4>i + bj-b 



A ' 



(3.17) 



where • is the sesquilinear dot product. We consider the right hand side of equation (3.17) as a 
function of A. As Bf = 0, this expression is a linear combination of the functions B^4)^ )^ of A, 
with coefficients given by the complementary boundary forms. 
The definitions of B^. and (f)\ imply 



BUt = 



1 



27rA(A) 

^ ' 1=1 3=1 



5;]^detX'^(A)M+(A)B,*(e— '"^^O 



n n 



E E X' ' {X)M+ {\)Mi u (A). 

^ ' 1=1 3=1 



But 



so 



5;]detX'^(A)Afifc(A) = A(A),5. 



3 



1 = 1 



BUI = ^^tkW- 



Similarly, 



B\ 



27r 



Finally, by equations (2.12), are polynomials of order at most n — 1. 



(3.18) 
(3.19a) 

(3.19b) 

□ 
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of Proposition 3.2 Let q be the solution of the problem. Then, since q satisfies the partial differ- 
ential equation (2.11a), 

-^^F+{q{;t}) = -aF+{S{q{;t))) = -aA"F+(g(.,t)) - aP+._,)(A), (3.20) 

where, by Lemma 3.3 is a polynomial of degree at most n — 1. Hence 



d 
dt 



Integrating with respect to t and applying the initial condition (2.11b), we find 



Similarly, 



F^{q{;t))^e'^^"'F^{f) 



(3.21) 



(3.22) 



(3.23) 



where P^^ is another polynomial of degree at most n — 1. The validity of the type II transform 
pair. Proposition 2.3 implies 



q{x,t) 



r+ 



P+(/)dA+ / e*^(^-i)-''-^"*P;^(/)dA 



IT- 
^iXx—aX^t 

^iX(x—l) — aX^t 

AXx—aX^t 



^iX{x~l)--aX^t 



j\^^"^p-,^^^iX)ds^ dX. (3.24) 



As are polynomials, the integrands 



iXx — aX^t 



^*e'^^"^P+ ^^)(A)ds^ and e^^(^-i)-'^^"* ' p-^, ^^{X) dt 



are both entire functions of A. Hence the third and fourth terms of equation (3.24) vanish. Inte- 
gration by parts yields 



AXx—aX^t 



^iX{x—l) — aX^t 



^ax^sp+ (A)ds ) = 0{X-') as A ^ oo 



e '^g(.,s)(-^)ds ) = 0{X-') as A -> oo 



within the region 
enclosed by r+, 

within the region 
enclosed by F". 



Hence, by Jordan's Lemma, the final two terms of equation (3.24) vanish. 



□ 



Remark 2. The same method may be used to solve initial-boundary value problems with inho- 
mogeneous boundary conditions. The primary difference is that statement (2.18) must be replaced 
with [ini Lemma 4.1]. 
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4 Analysis of the transform pair 



In this section we analyse the spectral properties of the transform pairs using the notion of aug- 
mented eigcnfunctions. 



4.1 Linearized KdV 
Augmented Eigenfunctions 

Let S"*^^-* and S^^^^ be the differential operators representing the spatial parts of the IBVPs 1 and 2, 
respectively. Each operator is a restriction of the same formal differential operator, {—id/ dx)^ to 
the domain of initial data compatible with the boundary conditions of the problem: 



= {/ e C^[0, 1] : /(O) = /(I) = 0, /'(O) = 2/'(l)}, 
2?(^(")) = {/ e C°-[0, 1] : /(O) = /(I) = /'(I) = 0}. 



(4.1) 
(4.2) 



A simple calculation reveals that {Fx : A G TQ,A{a) — 0} (where F\ is defined by equa- 



tions (1.8a), (1.8c) and (1.8dl) is a family of type I augmented eigenfunctions of S**^^^. Indeed, 



integration by parts yields 



i^A(5m/) 




(4.3) 



For any /, the remainder functional is an entire function of A and To is a closed, circular contour 



hence ( 1.12 ) holds 



In the same way {Fx : A G Fq} (where F\ is defined by equations (1.8al, (1.8c) and (1.8d)) is 
a family of type I augmented eigenfunctions of S^-^^\ Indeed 



X'Fxif) 
X'Fxif) 



f /"(O) - ^/'(O) 



2tt 



/"(I) 



A G 
A G 



(4.4) 



so the remainder functional is again entire. 

Furthermore, the ratio of the remainder functionals to the eigenvalue is a rational function with 



no pole in the regions enclosed by F^ and decaying as A — ^ oo. Jordan's lemma implies (1.13) 



hence {Fa : A G F+ U F } is a family of type II augmented eigenfunctions of the corresponding 
or 



Spectral representation of S'*^^^^ 

We have shown above that {Fx : A G Fq} is a family of type I augmented eigenfunctions and 
{Fx : A G F+ U F~} is a family of type II augmented eigenfunctions of S^^^\ each with eigenvalue 
A"^. It remains to show that the integrals 



e'^^Fx{Sf)dX, 



e'^^Fxif)dX 



(4.5) 



r+ur- 



converge. 

A simple calculation reveals that Fx{ip) has a removable singularity at A = 0, for any ip G C. 
Hence the first integral not only converges but evaluates to 0. Thus, the second integral represents 



fxiFxif)) = f and converges by Proposition 2.1 



This completes the proof of Theorem |1.5| for problem 2 
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Spectral representation of S^^'' 



By the above argument, it is clear that the transform pair (Fx, fx) d efined by equations (1.3) 
provides a spectral representation of S^^^ in the sense of Definition |l.4[ verifying Theorem 1.5 for 
problem 1. 

It is clear that {Fx : A G F^} is not a family of type I augmented eigenfunctions, so the 



1.3 



representation ( |1.3[ ) does not provide a spectral representation of 5*^^^ in the sense of Definition 
However, equation ( |1.7[ ) does provide a representation in the sense of Definition 1.3 Indeed, 
equation (1.7) implies that it is possible to deform the contours onto 



u 



creC: 
A(o-)=0 

It is possible to make this deformation without any reference to the initial-boundary value problem. 
By an argument similar to that in the proof of Proposition 2.1, we are able to 'close' (whereas in 
the earlier proof we 'opened') the contours F* onto simple circular contours each enclosing a single 
zero of A. Thus, an equivalent inverse transform is given by (1.9 1. It is clear that, for each a a 
zero of A, {F\ : A G Fo-} is a family of type I augmented eigenfunctions of S^^^ up to integration 
over Fcr. 

It remains to show that the series 



E 

crSC: 
A(cr)=0 



FxiSf) d\ 



(4.6) 



converges. The validity of the transform pair {F\, f^) defined by equations (1.8a), (1.8c), (1.8d) 



and ( 1.9 ) is insufficient to justify this convergence since, in general, S f may not satisfy the boundary 



conditions, so S f may not be a valid initial datum of the problem. Thus, we prove convergence 
directly. 

The augmented eigenfunctions F\ are meromorphic functions of A, represented in their defini- 



tion (1.8a), (1.8c), (1.8d) as the ratio of two entire functions, with singularities only at the zeros 



of the exponential polynomial A. The theory of exponential polynomials [7] implies that the only 
zeros of A are of finite order, so each integral in the series converges and is equal to the residue 
of the pole at a. Furthermore, an asymptotic calculation reveals that these zeros are at 0, a-'Afe, 
^k, for each j G {0, 1, 2} and fc G N, where 



A,. 



2A: - i j TT + i log 2 + O (e-^'^'^) 



— y/Skn 



fik = - [2k - -]7r + i\og2 + 0{e 



Evaluating the first derivative of A at these zeros, we find 

A'(Afe) = (-l)'=+iy2e''^'°s2gy3.(fc-i/6) ^O(^)^ 
A'(^fe) = (-l)'=\/2e'^^'°s2gy3.(fc-i/6) ^O(^)^ 



(4.7) 
(4.8) 



(4.9) 
(4.10) 



Hence, at most finitely many zeros of A are of order greater than 1. A straightforward calculation 
reveals that is a removable singularity. Hence, via a residue calculation and integration by parts, 
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we find that we can represent the tail of the series (4.61 in the form 

oo . 



k=N 



AfeA'(Afe) 



+ a^e'^^''- ((5/)(l)ri(aAfe) - (5/)(0)ro(«Afc)) 
[e^^-^ {{Sf){l)YM - (SfmYoipk)) 



((5/)(l)ri(aVfe) - (5/)(0)yo(aVfc))] + 0(fc"')} , (4.11a) 



where 





= 3 + 2(a2_^)g«aA 


+ 2(a-l)e*"'\ 


i"o(A) 






Zi(A) 


ae*"^ + 2e-*"^ + 




^o(A) 


6+ (a^ - l)e-*"^ 


+ (a - l)e-^"'^. 



(4.11b) 
(4.11c) 
(4.11d) 
(4.11e) 

As y, , e 0(exp(\/37rfc)), the Ricmann-Lebesgue lemma guarantees conditional convergence for 
all X e (0,1). 

This completes the proof of Theorem |1.6| 



Remark 3. We observed above that is removable singularity of F\ defined by (1.8a), (1.8c 



and (|l.8d|). The same holds for F\ defined by (1.8a), (1.8e) and (1.8f). Hence, for both problems 1 

e*^'^FA(/)dA = (4.12) 



and 2, 



and we could redefine the inverse transform (|1.8b|) as 
^(A) ^ /(x) : 



y +J |e'^"i^(A)dA, xe [0,1] 



(4.13) 



This permits spectral representations of both S^^^ and S^^^^ via augm ented eigenfunctions of type II 
only, that is spectral representations in the sense of Definition 1.4 but with iS^^' = 0. 

4.2 General 



We will show that the transform pair {Fx, fx) defined by equations (2.14) represents spectral 
decomposition into type I and type II augmented eigenfunctions. 

Theorem 4.1. Let S be the spatial differential operator associated with a type II IBVP. Then the 
transform pair {Fx, f^) provides a spectral representation of S in the sense of Definition 1 



The principal tools for constructing families of augmented eigenfunctions are Lemma |3.3[ as 
well as the following lemma: 

Lemma 4.2. Let F^ be the functionals defined in equations ( |2.15 ). 

(i) Let 7 be any simple closed contour. Then {F^ : A G 7} are families of type I augmented 
eigenfunctions of S up to integration along 7 with eigenvalues A" . 
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(ii) Let 7 be any simple closed contour which neither passes through nor encloses 0. Then {F^ : 
A € 7} are families of type II augmented eigenf unctions of S up to integration along 7 with 
eigenvalues A". 

(Hi) Let ^ < 0' ^ TT and define 7^ to be the boundary of the open set 

{A e C : |A| > e, 61 < argA < 6*'}; (4.14) 

similarly, 7" is the boundary of the open set 

{A e C : |A| > e, -6*' < argA < -6*}. (4.15) 

Both 7+ and 7^ have positive orientation. Then {F^ : A G 7*} are families of type II 
augmented eigenf unctions of S up to integration along 7* with eigenvalues A" . 

Proof. 



(i) & (ii) By Lemma 3.3 the remainder functionals are analytic in A within the region bounded 
by 7. Cauchy's theorem yields the result. 



(iii) The set enclosed by 7^ is contained within the upper half-plane. By Lemma 



3.3 



J'Y+ 



A-"(F+(5/)-A"^+(/))dA= / e'^-A-"P/(A)dA, (4.16) 



and the integrand is the product of e*^^ with a function analytic on the enclosed set and 
decaying as A — ?> 00. Hence, by Jordan's Lemma, the integral of the remainder functionals 
vanishes for all x > 0. For 7", the proof is similar. 

□ 

Remark 4. If we restrict to the case < 9 < 9' < tt then the functionals F^ form families of type I 
augmented eigenfunctions up to integration along the resulting contours but this is insufficient for 
our purposes. Indeed, an infinite component of Fq lies on the real axis, but 

00 

e*^^P+(A)dA (4.17) 

-co 

diverges and can only be interpreted as a sum of delta functions and their derivatives. 

Let (5*, a) be such that the associated initial-boundary value problem is well-posed. Then there 
exists a complete system of augmented eigenfunctions associated with S, some of which are type I 
whereas the rest are type II. Indeed: 

Proposition 4.3. The system 

-Fo = {F+ : AeF+}U{^^- : AeF„-} (4.18) 

is a family of type I augmented eigenfunctions of S up to integration overTg, with eigenvalues A". 
The system 

Ta = {F+ :XeT+}U{F^ -.XeT-} (4.19) 

is a family of type II augmented eigenfunctions of S up to integration over Fq, with eigenvalues 
A". 

Furthermore, if an initial-boundary value problem associated with S is well-posed, then T = 
J^o U J^Q is a complete system. 
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Proof. Considering / e <I> as the initial datum of the homogeneous initial-boundary value problem 
and applying Proposition 3.2 we evaluate the solution of problem ( |2.11 1 at t = 0, 



dX. 



(4.20) 



Thus, if F± (/) = for all A G Tq then / = 0. 



By Lemma 4.2 (i), J^a is a system of type I augmented eigenfmictions up to integration along 



r+uTo". 



Applying Lemma 3.3 to J-a, we obtain 

A 



F±(5/) = A"F±(/) + i?±(/), 



with 



(4.21) 



(4.22) 



4.2 



(ii), we can deform the contours onto the union of several contours of the form 



By Lemma 

of the 7^ appearing in Lemma 4.2 (iii). The latter result completes the proof. 



of Theorem ^.1 



equations (1.19 



□ 



Proposition |4.3| establishes completeness of the augmented eigenfunctions and 
, under the assumption that the integrals converge. The series of residues 



/ e'^''F^{Sf)d\ = 2TTi y e^-^^Resi^, 

/r' ' ^ A— cr 



'■{Sf). 



(4.23) 



A(<t)=0 



whose convergence is guaranteed by the well-posedness of the initial-boundary value problem [TT] . 
Indeed, a necessary condition for well-posedness is the convergence of this series for Sf S <&. But 
then the definition of implies 

ResF±(/) = 0{\a\-'-^), where j = max{A: : V / G $, /(^^(O) = /C^^l) = 0}, 

A — (7 

so ResA=cr ^a(<S'/) = 0(|cr|^^) and the Riemann-Lebesgue lemma gives convergence. This verifies 
statement (1.17). Theorem 2.3 ensures convergence of the right hand side of equation (1.19b|. 
Hence statement ( |1.18 1 holds. □ 

Remark 5. Suppose S* is a type I operator. 

By the definition of a type I operator (more precisely, by the properties of an associated type I 
IBVP, see [U), F^^cj)) = 0(A~^) as A — )■ oo within the sectors interior to Pj. Hence, by Jordan's 
Lemma, 



^iXx 77' + 



^A^WdA 



^-i^-(0)dA = O. 



Hence, it is possible to define an alternative inverse transform 
^^(A) ^ fix) : f^{F) = 



F{X) dA, 



(4.24) 



(4.25) 



equivalent to fx- The new transform pair (Fx, /J") defined by equations (2.14) and (4.25) may be 
used to solve an IBVP associated with S hence 



(4.26) 



is a complete system of functionals on $. 

Moreover, J^) is a family of type I augmented eigenfunctions only. Hence, J^o provides a spectral 
representation of S in the sense of Definition [L3] Via a residue calculation at each zero of A, one 
obtains a classical spectral representation of 5 as a series of (generalised) eigenfunctions. 

We emphasize that this spectral representation without type II augmented eigenfunctions is 
only possible for a type I operator. 
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Remark 6. By definition, the point 3e/2 is always exterior to the set enclosed by F. Therefore 
introducing a pole at Zejl does not affect the convergence of the contour integral along F. This 
means that, the system T' = {(A — 3e/2)~"Fx : A G F} is a family of type I augmented eigenfunc- 
tions, thus no type II augmented eigenfunctions are required; equation (1.16) holds for T' and the 



integrals converge. However, we cannot show that T' is complete, so we do not have a spectral 
representation of S through the system JF' . 

Remark 7. There may be at infinitely many circular components of Fq, each corresponding to 
a zero of A which lies in the interior of a sector enclosed by the main component of F^. It is 



clear that in equations (2.17) and (3.12), representing the validity of the transform pair and the 



solution of the initial-boundary value problem, the contributions of the integrals around these 
circular contours are cancelled by the contributions of the integrals around certain components 
of Fq, as shown in Figure [3] Hence, we could redefine the contours F^ and Fo to exclude these 
circular components without affecting the validity of Propositions |2.3| and |3.2[ 

Our choice of F^ is intended to reinforce the notion that S is split into two parts by the 
augmented eigenfunctions. In Fq, we have chosen a contour which encloses each zero of the charac- 
terstic determinant individually, since each of these zeros is a classical eigenvalue, so J^o corresponds 
to the set of all generalised eigenfunctions. Hence Ta corresponds only to the additional spectral 
objects necessary to form a complete system. 

Remark 8. As Fa encloses no zeros of A, we could choose a i? > and redefine F^^ as the 
boundary of 

{A e : |A| > i?, Re(aA") > O} \ |J D{a, 2e), (4.27) 

o-gC: 
A(cr)=0 

deforming F^ over a finite region. By considering the limit R — > oo, we claim that J-^ can be seen 
to represent spectral objects with eigenvalue at infinity. 



Remark 9. By Lemma 



4.2 



ii), for all cr 7^ such that A{a) = 0, it holds that {F^ : A G C(cr, e)} 



are families of type II augmented eigenfunctions. Hence, the only component of Fq that may not 
be a family of type II augmented eigenfunctions is C{0,e). If 



crGC+: 
A((t)=0 



7a"-r-U y C{a,e), 
ere 

Ai 

7o = C(0,£), 



A{a)=0 



then 



J-: = {i^+A:Ae7+}U{i^-A: Ae7-} 
is a family of type II augmented eigenfunctions and 

J-^ = {F+X : A e 7o} 



(4.28a) 

(4.28b) 
(4.28c) 
(4.29) 

(4.30) 



is a family of type I augmented eigenfunctions of S. 



a spectral representation of S in the sense of Definition 1.4 



For S type I or type II, F'^ U Fq provides 
with minimal type I augmented 



eigenfunctions. (Note that it is possible to cancel certain circular components of 7^ •) 

Assume that is a remo vabl e singularity of F^ . Then F'^ provides a spectral representation of 



S in the sense of Definition 1.4 with E^^^ = 0. We have already identified the operators S'--^^ and 



5'(ii) for which this representation is possible (see Remark [3]) 
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Remark 10. The vaHdity of Lemmata 3.3 and 4.2 does not depend upon the class to which S 
belongs. Hence, even if all IBVPs associated with S are ill-posed, it is still possible to construct 
families of augmented eigenfunctions of S. However, without the well-posedness of an associated 
initial-boundary value problem, an alternative method is required in order to analyse the complete- 
ness of these families. Without completeness results, it is impossible to discuss the diagonalisation 
by augmented eigenfunctions. 



5 Conclusion 

In the classical separation of variables, one makes a particular assumption on the form of the 
solution. For evolution PDEs in one dimension, this is usually expressed as 

"Assume the solution takes the form q{x,t) = T{t)^{x) for all {x,t) € [0,1] x [0,r] for 
some C G C°°[0, 1] and r € C°°[0,T]." 

However, when applying the boundary conditions, one superimposes infinitely many such solutions. 
So it would be more accurate to use the assumption 

"Assume the solution takes the form q{x,t) = J2m<£N'^m{t)^m{x) for some sequences of 
functions € C°°[0, 1] which are eigenfunctions of the spatial differential operator, 
and r„i e C°° [0, T] ; assume that the series converges uniformly for {x,t) € [0, 1] x [0, T] ." 

For this 'separation of variables' scheme to yield a result, we require completeness of the eigen- 
functions {^m)meN in the space of admissible initial data. 

The concept of generalized eigenfunctions, as presented by Gelfand and coauthors |5l |6j allows 
one to weaken the above assumption in two ways: first, it allows the index set to be uncountable, 
hence the series is replaced by an integral. Second, certain additional spectral functions, which are 
not genuine eigenfunctions, are admitted to be part of the series. 

An integral expansion in generalized eigenfunctions is insufficient to describe the solutions of 
IBVPs obtained via the unified transform method for type II problems. In order to describe these 
IBVPs, we have introduced type II augmented eigenfunctions. Using these new eigenfunctions, the 
assumption is weakened further: 

"Assume the solution takes the form q{x, t) — /^^^p T„i{t)^„i{x) dm for some functions 
£,m G C°°[0, 1], which are type I and II augmented eigenfunctions of the spatial differ- 
ential operator, and S C°°[0,T]; assume that the integral converges uniformly for 

{x,t) e [0,1] X [o,r]." 

It appears that it is not possible to weaken the above assumption any further. Indeed, it has 
been established in |4] that the unified method provides the solution of all well-posed problems. 
The main contribution of this paper is to replace the above assumption with the following theorem: 

"Suppose q{x, t) is the C°° solution of a well-posed two-point linear constant-coefficient 
initial-boundary value problem. Then q{x,t) = J^^p Tm(i)Cm(a^) dm, where £ 
C°°[0, 1] are type I and II augmented eigenfunctions of the spatial differential operator 
and T,„ G C°°[0,T] are some coefficient functions. The integral converges uniformly for 

{x,t) e [0,1] X [o,r]." 



In summary, both type I and type II IBVPs admit integral representations like ( 1.3 1, which give 
rise to transform pairs associated with a combination of type I and type II augmented eigenfunc- 
tions. For type I IBVPs, it is possible (by appropriate contour deformations) to obtain alternative 



integral representations like (1.7), which give rise to transform pairs associated with only type I 
augmented eigenfunctions. Furthermore, in this case, a residue calculation yields a classical series 
representation, which can be associated with Gel'fand's generalised eigenfunctions. 

The research leading to these results has received funding from the European Union's Seventh 
Framework Programme (FP7-REGPOT-2009-1) under grant agreement n° 245749. 
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